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We employ the dynamical mean field approximation to perform a systematic study of magnetism 
in Gai_a;Mni:As. Our model incorporates the effects of the strong spin-orbit coupling on the J = 3/2 
GaAs valence bands and of the exchange interaction between the randomly distributed magnetic 
ions and the itinerant holes. The ferromagnetic phase transition temperature Tc is obtained for 
different values of the impurity-hole coupling Jc and of the hole concentration nn at the Mn doping 
of a: = 0.05. We also investigate the temperature dependence of the local magnetization and spin 
polarization of the holes. By comparing our results with those for a single band Hamiltonian in which 
the spin-orbit coupling is switched off, we conclude that the spin-orbit coupling in Gai-ajMn^^As 
gives rise to frustration in the ferromagnetic order, strengthening recent findings by Zarand and 
Janko (Phys. Rev. Lett. 89, 047201 (2002)). 



The discovery of ferromagnetism in GaAs doped with 
Mn has renewed interest in the properties of diluted mag- 
netic semiconductors. Q Since these materials are good 
sources of polarized holes, they may form the basis of a 
new field called spintronics. 2] Spintronic devices employ 
both the spin and the charge of the carrier to convey 
information. 

In Gai_2,Mna;As, the Mn ions are in the Mn^+ state 
with a half-filled d shell of total spin S = 5/2. HH Since 
Mn^"*" ions primarily replace Ga'^"'", they act as effective 
acceptors by supplying holes as well as localized spins. 
The valence band of pure GaAs is p-like so the strong 
spin-orbit interaction couples the I = 1 angular momen- 
tum of the p orbitals to the electron spin (s = 1/2), 
resulting in a total spin 7 = ^-1-5 = 3/2 for the valence 
holes. 15| As discussed by Zarand and Janko, [cj the strong 
spin-orbit coupling also induces an anisotropic carrier- 
mediated interaction between the Mn ions and, as a con- 
sequence, frustration in their ferromagnetic order. While 
the results of Ref. were limited to the metallic regime 
and to small values of the Mn-hole coupling, a different 
approach suggests the presence of an impurity band in 
the dilute limit. 

In this Letter, we employ the Dyiiamical Mean Field 
Approximation (DMFA) 8,, H US [iJl to perform a sys- 
tematic analysis of ferromagnetism in Gai_a;Mna;As, in- 
cluding the effects of strong spin-orbit coupling on the 
J = 3/2 GaAs valence band. The DMFA includes the 
spin-split impurity band through quantum self-energy 
corrections which are not included in other mean-field 
theories. Because this method is non-perturbative, it al- 
lows us to study both the metallic and impurity-band 
regimes as well as both small and large couplings. We 
show how the spin-orbit interaction affects the ferro- 
magnetic critical temperature Tc, the hole polarization, 
and the Mn magnetization. By comparing our results 
with those for a single band Hamiltonian without spin- 
orbit coupling, we conclude that strong spin-orbit cou- 



pling in Gai_a:Mna;As produces frustration for all cou- 
pling strengths. For carrier concentrations smaller than 
the doping, both Tc and the polarization of the carriers 
are reduced for all values of the coupling. For larger car- 
rier densities and large couplings, we unexpectedly find 
that frustration induces a small but finite Tc, in sharp 
contrast with the vanishing Tc found when the spin-orbit 
is neglected. 

Our starting point is the simplified Hamiltonian pro- 
posed in Ref. 

H = Ho-JcY.S,-J{R,). (1) 

The first term includes the electronic dispersion and the 
spin-orbit coupling of the J ^3/2 valence holes within 
the spherical approximation. 12] The second term is the 
dominant part of the interaction between the Mn spins 
and the valence holes, with Jc the exchange coupling 
and 3{Ri) the total J = 3/2 spin density of the holes at 
the site i of a Mn ion with spin S^. The relatively large 
magnitude of the Mn local moment justifies a classical 
treatment of its spin. 

Within the spherical approximation, the Hamiltonian 
of pure GaAs is rotationally invariant. Hence, Hq is di- 
agonal in a chiral basis, 

^o = E^4.5.., (2) 

where cj^ ^ creates a chiral hole with momentum k paral- 
lel to its spin and J-k = ±3/2 or ±1/2. The two effective 
masses = w 0.5to and — mi k, O.OTto corre- 
spond to the heavy and light bands with 7 = ±3/2 and 
±1/2 respectively (m is the electron mass). 

Unfortunately, the exchange interaction term obtains 
a rather complicated momentum-dependent form in the 
chiral basisfq that is responsible for the frustrated order 
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of the Mn. The competition between the strong spin- 
orbit couphng on the hole bands, which ahgns the hole 
spin parallel to its momentum, and the exchange inter- 
action with the local moment, which aligns the hole and 
local spins, precludes all of the carrier density from me- 
diating the magnetic order. 

We develop a DMFA algorithm that takes advantage of 
the simple diagonal form of Hq in the chiral basis and the 
local form of the exchange interaction in the non-chiral 
basis. The coarse-grained Green function matrix in the 
non-chiral fermion basis is 



G{iLLln) 
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(3) 



where N is the number of k points in the first Bril- 

louin zone and e(fc) = R^k) R{k) is the disper- 

2m-y 

sion in the spherical approximation. Here, R are spin 
3/2 rotation matrices that relate the fermion operator, 
to its chiral counterpart, Ck'y = R-y^ik^Ck^. The 
mean-field function Q{iujn) = [G^^iitOn) + S(ict;„)]~^ is 
required to solve the DMFA impurity problem. At a 
nonmagnetic site, the local Green function equals the 
mean-field function, Gnon = Q- At a magnetic site, Gmg 
must solve a local problem. By treating disorder in a 
fashion similar to the coherent potential approximation 
(CPA) 14] for a given local spin configuration, we obtain 

Now Gmg{ioJn) must be averaged over all possible spin 
orientations at the local site and over all possible impu- 
rity configurations on the lattice. The former is imple- 
mented by using the effective action 



Seff{S) = - ^logdet \g{iUr,){g'\lUJn) + JcS • J) 



xe*'^"° (4) 
to average over the angular distribution of the local spins: 

GmgiiuJn)'^ = ^ J dfls G„ig(iw„)exp[-S'e//(S)] , 

(5) 

where Z = / df^s exp[-S'e//(S)]. If the Mn 
ions are randomly distributed with probability x, 
the configurationally-averaged Green function reads 

Gavg{i(^n) = (Gmg{itOn)^ X + ^/(za;„)(l - x) . 

Finally, the magnetization of the Mn ions can be cal- 
culated as: 



(6) 



where a small magnetic field SH^ is applied to break 
the symmetry along a preferential direction, i.e., the z 
axis, [iai By fitting the magnetization (T) in the vicin- 
ity of the transition to the Curie- Weiss form, we can ex- 
tract the value of Tc for each set of parameters studied. 
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FIG. 1: Phase transition temperature Tc versus carrier con- 
centration Uh for Mn doping x = 0.05, various Jc values and 
effective mass ratios, a. 



We focus on the doping of x — 0.05 associated with 
the highest reported. [iH El IHliS 113 In order to 
clearly elucidate the role played by the spin-orbit cou- 
pling, we introduce two parameters. One is the hole 
mass of an equivalent system composed of two degen- 
erate bands with the same average kinetic energy as our 

system: 2miq = nih ^ + mi^ . The other is the ratio of 

Tni 

light and heavy hole effective masses, a = . The chi- 

rality can be switched off by setting a — 1 while keeping 
nieq constant. 

Fig.^presents our results for the phase transition tem- 
perature Tc as a function of the hole concentration Uh for 
a fixed Mn concentration, x — 0.05, and for various val- 
ues of Jc and a. There are two regimes corresponding to 
small and large values of the interaction strength, Jc- For 
large Jc, the holes form bound states with the Mn impu- 
rities, an impurity band develops inside the GaAs gap, 
and the properties of the host are greatly afi^ected. The 
value of Jc at which the impurity band appears depends 
upon the value of a. 

First consider the non-chiral case with a = 1.0 so that 
the spin-orbit coupling is turned off. For Jc = 1.0 eV 
the system is far from the impurity-band regime and Tc 
has a relatively slow variation with respect to rih- For 
Jc = 2.0 eV the system is beyond the threshold for the 
formation of an impurity band, which then dominates 
the physics. The maximum Tc occurs when the impurity 
band is nearly half filled {uh « x/2).^ For Jc = 4.0 
eV the impurity band is well established and the maxi- 
mum Tc is large. Because ferromagnetic order restricts 
the hopping of holes when the impurity band is full, Tc 
vanishes for Jc > 1.0 eV and Uh > x. An antiferromag- 
netic ground state is energetically more favorable in the 
regime Uh > x because the carriers can then easily 
hop from one impurity site to another. Magnetization 
curves for J = 2.0 and 4.0 eV and for rih > x provide 
evidence for antiferromagnetism. They fit a Curie- Weiss 
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FIG. 2: Chemical potential, ^, versus Uh for Mn doping x = 
0.05 and T = 0.05 eV using the same values of Jc and a as in 
Fig. 1. 



functional form, A'P = —, — with a positive O rane- 

l + T/e ^ 

ing from 40 to 117K, indicating antiferromagnetic order 
for a = 1.0 and Uh > x. 

We now consider the chiral case with a = 0.14 so that 
the spin-orbit coupling in GaAs is turned on. The effects 
of frustration on the ferromagnetic critical temperature 
are readily seen in Fig. For Jc = 1.0 eV, Tc consis- 
tently lies below its non-chiral counterpart. At Jc = 2.0 
and 4.0 eV, the difference between the chiral and non- 
chiral results is even larger. For < x, the non-chiral 
Tc is always higher than the chiral Tc- However, for 
nil > X, the chiral results continue to yield a finite al- 
though greatly diminished Tc in the same regime where 
Tc vanished when a — 1.0. Surprisingly, the ferromag- 
netic Tc survives at large dopings due to the intrinsic 
frustration in the system. Strong spin-orbit scattering 
allows the impurity-band carriers to hop between impu- 
rity sites that are ferromagnetically aligned even when 
the impurity band is full, thereby stabilizing the mag- 
netic order. 

The formation of an impurity band is responsible for 
many of the outstanding properties of this svstem.p3| 
Impurity-band formation appears as a discontinuity of 
the hole chemical potential fi at Uh = x. Fig. [21 depicts /i 
versus Uf^ for T = 0.05 eV {580K) and for the same values 
of Jc and a used in Fig. ^ The presence of an impurity 
band depends on both Jc and a. When Jc — 1.0, no 
impurity band is present so the chiral and non-chiral fi 
are almost identical and neither shows a discontinuity. 
When Jc = 2.0 eV and a = 1.0, the impurity band has 
already split from the valence band. At the same value of 
Jc — 2.0 eV but when a = 0.14, the change in slope of /i 
aX Uh — X indicates that the impurity band is present but 
that it overlaps with the main band. When Jc — 4.0 eV 
and a = 0.14, the hole chemical potential is discontinuous 
at nh = X, signaling the splitting of the impurity band 
from the valence band. For these larger values of Jc, the 
non-chiral always lies below the chiral value, indicating 
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FIG. 3; Magnetization versus temperature for different chi- 
rality modes and different values of Jc at = x/2 = 0.025 
and 5H^ — 0.004 eV. The curves are the static mean field fits 
with Tc values equal to those given by our results, solid curves 
fit the results for a = 0.14; dashed curves the ones for a = 1. 



that the non-chiral impurity band lies at lower energies 
than the chiral one. Since the interaction between the 
impurity moments is mediated by the host of holes, the 
less pronounced impurity band for a = 0.14 is another 
signature of the frustration produced by chirality. 

Optical conductivity measurements j23| show that the 
impurity band is split from the valence band at the dop- 
ing of a; = 0.05. This fact and our results for Tc suggests 
a value of Jc between the value of 1.2 eV obtained from 
photoemission techniques^ and the value of 4.5 ~ 6 eV 
inferred from infrared spectroscopy. |3 

The magnetization also reveals the effects of frustra- 
tion. In Fig. the temperature dependence of the mag- 
netization is plotted with Uh — x/2 — 0.025. Our results 
are compared with Curie- Weiss fits for the same Tc val- 
ues. When Jc = 1.0 eV, the mean-field curves perfectly 
fit the data not only at large temperatures but also below 
Tc- But for larger Jc values and low temperatures, dis- 
crepancies appear between the static mean-field curves 
and the DMFA results. For a = 1.0, the magnetization 
in our model lies above the mean-field Heisenberg mag- 
netization when T < Tc since the scattering between the 
itinerant carriers and the localized spins becomes coher- 
ent at low temperatures, thereby enhancing the magnetic 
order of the local ions. This effect has also been seen in 
the double exchange model. 24] However, for a = 0.14, 
frustration reduces the low-temperature coherence be- 
tween the carriers and local ions. So when T < Tc, the 
magnetization curves lie below the mean-field Heisenberg 
predictions. 

Lastly, we study the carrier spin polarization, V, which 
is constructed from appropriate sums of the Green func- 
tion. Control of V is important for spintronics appli- 
cations, since the spin polarization of the holes is re- 
quired to transport information. Fig. ^ illustrates the 
temperature dependence of the hole polarization when 
n-h = x/2 = 0.025 for several values of Jc and a. For 
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FIG. 4: The hole polarization P versus temperature at rih — 
x/2 = 0.025 and 5H^ = 0.004 eV. 

all Jc values, chirality suppresses the polarization. For 
a ~ 1.0, large J^, and low temperatures, V approaches 
3/2 since all holes occupy the lowest energy level with 
= +3/2. But due to frustration even for large Jc 
and low temperatures, V for a = 0.14 is significantly 
smaller. This agrees with previous calculations of the 
zero-temperature polarization, where the destructive ef- 
fects of the spin-orbit interaction were also found. 0| 

In summary, we have studied the diluted magnetic 
semiconductor Gai^^Mn^As using a non-perturbative 
multi-orbital DMFA algorithm incorporating the effects 
of the strong spin-orbit coupling. We have calculated the 
ferromagnetic critical temperature, hole chemical poten- 
tial, local ion magnetization, and hole polarization for a 



broad range of model parameters. We find that the spin- 
orbit coupling leads to frustration and reduced magne- 
tization when the hole concentration Uh is smaller than 
the impurity concentration x, in agreement with previ- 
ous perturbative calculations. 6] In addition, we find that 
this behavior persists for large values of Jc, and that frus- 
tration greatly reduces the transition temperature Tc and 
the polarization of the carriers for all Jc- Finally, when 
Jc is large, we find the surprising result that frustration 
induces a region of finite Tc for rih > x. 

This approach has promising future ramifications. It 
can be extended to study other magnetic semiconductors 
and realistic devices such as semiconducting heterostruc- 
tures and quantum dot systems, which can be tailored to 
take full advantage of the intrinsic anisotropy of the fer- 
romagnetic order. More sophisticated approaches, such 
as the dynamical cluster approximation (DCA).)2fi] may 
be used to go beyond the single-site approximation and 
explore the cooperative and glassy effects of frustration, 
such as the reduction in the local magnetization at low 
temperatures. 
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